Band-structure topologies of graphene: spin-orbit coupling effects from first principles 
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The electronic band structure of graphene in the presence of spin-orbit coupling and transverse 
electric field is investigated from first principles using the linearized augmented plane-wave method. 
The spin-orbit coupling opens a gap at the K{K')-pomt of the magnitude of 24 fieV (0.28 K). 
This intrinsic splitting comes 96% from the usually neglected d and higher orbitals. The electric 
field induces an additional (extrinsic) Bychkov-Rashba-type splitting of 10 /xeV (0.11 K) per V/nm, 
coming from the a-n mixing. A 'mini-ripple' configuration with every other atom is shifted out of 
the sheet by less than 1% differs little from the intrinsic case. 
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The fascination with graphene [l[ , the one-atom-thick 
allotrope of carbon, comes from its two-dimensional 
structure as well as from its unique electronic properties 
0, S 0, 0, 0| • The latter originate from the specific elec- 
tronic band structure at the Fermi level: electrons move 
with a constant velocity, apparently without mass and 
a spectral gap. Analogy with massless Dirac fermions 
is often drawn, presenting graphene as a solid-state toy 
for relativistic quantum mechanics. Ironically, this nice 
analogy is broken by the relativistic effects themselves. In 
particular, the interaction of the orbital and spin degrees 
of freedom, spin-orbit coupling, gives the electrons in 
graphene a finite mass and induces a gap in the spectrum. 
How large is the gap and which orbital states contribute 
to it? This question is crucial for knowing graphene's 
band-structure topology, understanding its spin trans- 
port and spin relaxation properties 0, 0| , or for assessing 
prospects of graphene for spin-based quantum computing 
(9|. By performing comprehensive first-principles calcu- 
lations we predict the spectral gap and establish the rel- 
evant electronic spectrum of graphene in the presence of 
external transverse electric field. We find that realistic 
electric fields can tune among different band structure 
topologies with important ramifications for the physics 
of graphene. 

Carbon atoms in graphene are arranged in a honey- 
comb lattice which comprises two triangular Bravais lat- 
tices; the unit cell has two atoms. The corresponding 
reciprocal lattice is again honeycomb, with two nonequiv- 
alent vertices K and K' which are the Fermi momenta 
of a neutral graphene. The states relevant for transport 
are concentrated in two touching cones with the tips at 
K{K') — the Dirac points — as illustrated in Fig. 1. The 
corresponding Bloch states are formed mainly by the car- 
bon valence Pz orbitals (the z-axis is perpendicular to the 
graphene plane) forming the two tt bands (cones). The 
other three occupied valence states of carbon form the 
deep-lying a bands by sp^ hybridization; these are re- 
sponsible for the robustness of graphene's structure. The 
states in the lower cones are hole or valence like, the up- 
per cone states are electron or conduction like, borrowing 
from semiconductor terminology. These essentials of the 



electronic band structure of grap hene were worked out 
many decades ago [13, El, Elil [3, [H ■ 

The above picture breaks down when spin-orbit cou- 
pling is included. The most important modification to 
the band structure is the opening of a gap at K{K'), 
as predicted by Slonczcwski and Weiss [IJ, [H. More 
severe changes occur when graphene is subject to trans- 
verse electric field that can come from the substrate or 
electric gates. Several questions arise: (i) What is the 
gap caused by the spin-orbit coupling itself? The mag- 
nitude of this intrinsic gap gives the temperature scale 
for observing the spin-orbit effects. While the spin-orbit 
coupling splitting of the p-orbitals in carbon atoms is of 
order 10 meV, the relevant states in graphene comprise 
Pz orbitals which have no net orbital momentum along z 
so the effect is expected to be rather weak, (ii) For real- 
istic electric fields, what are the corresponding extrinsic 
effects on the band structure? Are the effects compa- 
rable to the intrinsic ones? The answer helps to decide 
on the spin transport mode as well as on spin relaxation 
and spin coherence mechanisms, or to see if the band 
structure can be tailored, (iii) Which orbitals are in- 
volved in the intrinsic and extrinsic effects? Finally, (iv) 
do transverse deformations of the graphene plane induce 
significant spin-orbit effects? 

To investigate the electronic band structure of 
graphene in the presence of spin-orbit coupling and to 
answer questions (i) to (iv) we employed the full poten- 
tial linearized augmented plane waves (LAPW) method 
based on density functional theory [l3|. For exchange- 
correlation effects we used the generalized gradient ap- 
proximation (GGA) [l^. In our three-dimensional cal- 
culation the graphene sheets of lattice constant a = 
1.42-\/3A are separated by the distance 20 A, large 
enough for the inter-sheet tunneling to be negligible. In- 
tegration in the reciprocal space was performed by the 
modified Blochl tetrahedron scheme, taking the mesh of 
33 X 33 fc-points in the irreducible Brillouin zone wedge. 
As the plane-wave cut-off we took 9.87 A"-"^. The Is core 
states were obtained by solving the Dirac equation, while 
spin-orbit coupling for the valence electrons was treated 
within the muffin-tin radius of 1.34 a.u. by the second 
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variational method jl9| . Finally, external transverse elec- 
tric field was included as a periodic zigzag electric poten- 
tial added to the exchange-correlation functional [20[ • 

Our main results are shown in Fig. [2] which displays 
a variety of the band structure topologies of graphene, 
tunable by electric field. The intrinsic case (zero ex- 
ternal field) shows a splitting of the Dirac cones at the 
K{K') points into two "rounded" cones, with the gap of 
2A7 = 24 /^eV (0.28 K). The bands are two-fold degener- 
ate due to the presence of time reversal and space inver- 
sion symmetries pH . Transverse electric field breaks the 
latter symmetry, resulting in a spin-splitting 2Abr of the 
energy levels: for each momentum at a given band there 
are two states with energy differing by 2Abr. This ex- 
trinsic splitting is akin to the Bychkov-Rashba (B R) sp in- 
orbit coupling in semiconductor heterostructures |22l |. 

The band structure of graphene in the presence of 
transverse electric field depends rather strongly on the 
interplay of the intrinsic and extrinsic spin-orbit coupling 
effects. As the magnitude of the electric field increases, 
we encounter the topologies on display in Fig. O If the 
BR splitting is lower than the intrinsic one, the spectral 
gap gets smaller. The electron branch of the spectrum is 
still degenerate at K{K')\ in contrast, the hole branch 
is split by 2Abr. This topology can give a quantum 
spin Hall insulator Curiously, as the electric field 
is such that Abr ~ Ai, one of the hole branches rises, 
forming a genuine touching cones structure of massless 
fermions with one of the electron branches. The two 
remaining branches are parabolic (massive). For fields 
such that Abr > Ai, all the branches are again parabolic, 
with a degeneracy of one electron and one hole band. 
The calculated spectrum at K(K') follows the recipe 
M'^BR + J^IAbr — mAiI, with fj, and v being ±1. 

The physics behind the calculated spectral topolo- 
gies can be described qualitatively by previously pro- 
posed effective Hamiltonians. Without spin-orbit cou- 
pling the electronic band structure of graphene in the 
vicinity of K(K') is described by the Hamiltonian Ho = 
hv-p{K.'Jxkx + (Jyky). Here vp is the Fermi velocity, and 
ky are the cartesian components of the electron wave vec- 
tor measured from K{K'), the parameter k = 1 (—1) for 
the cones at K {K'), and ax and ay are the Pauli matri- 
ces acting on the so called pseudospin space formed by 
the two triangular sublattices of graphene. The Hamilto- 
nian Tio describes gapless states with conical dispersion 
Sq — vkvplkl near the Dirac points. The eigenstates are 
-0^ = I / ^/2\ve~^"-'^ , 1) for the electron band, v = \, 
and hole band, v = —1. 

The intrinsic spin-orbit coupling is described by the ef- 
fective Hamiltonian Hso = AiKtTzSz Heresz is the 
spin Pauli matrix. The spin-orbit coupling lifts the or- 
bital degeneracy at K{K'). Indeed, the eigenvalues of the 
combined Hamiltonian, Tio + Hso, are Si, = v\J eg -I- Aj . 
The bands are split by 2|Ai|, but the two- fold degen- 
eracy of the bands, required by space inversion and 



time reversal symmetry, remains. The eigenvectors are 
V'Ati^ = XAt|e"""^(K£i. + A«^i)/eo, ^)|C^,v■ Here X/x is the 
spin spinor with /i = ±1; the normalization constant is 
C^, = [l + (K£. + AiAi)Ve2]i/2_ 

The extrinsic spin-orbit coupling of the Bychkov- 
Rashba type in graphene can be described by the Hamil- 
tonian T^BR = ABR(tcra;Sy — aySx), where Abr is the 
Bychkov-Rashba parameter Q. Unlike in semiconduc- 
tor heterostructures, the coupling in graphene does not 
depend on the magnitude of the electron momentum, as 
the electrons at K{K') have a constant velocity. The 
electronic bands near K {K') are now modified to 



M^BR + i^V {hv^kY + (Abr - 



(1) 



The corresponding eigenvectors are 



= (x_|Ke-*"^[(e^,-Ai)/eo]", 1)+ (2) 
A^x+I - *A^e-'(i+-)'^, ie-'^[{h - £^..)/eo]"))/C^. , 

with the normalization constant C^^, = •\/2(l + [(Aj — 
^/ji>)/£o]^'^)^^- The expectation value of the spin. 
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\J^l + (^i - M^br)^ 




(3) 



is fc-dependent and lies in the graphene plane. The inclu- 
sion of the extrinsic coupling lifts the two-fold degeneracy 
of the bands. Only the time-reversal Kramers degeneracy 
remains, coupling states at K and K' . 

Our first-principles results show that the above effec- 
tive Hamiltonian model gives a remarkably faithful de- 
scription of graphenc's band structure at K{K'). The 
comparison is shown in Fig. [2] The dispersions given by 
Eq. Il]) differ from the numerical results by less than 5% 
up to ±200 meV away from the Fermi level. With the 
parameters supplied by the first-principles calculations 
the analytical model becomes highly accurate. 

The extrinsic splitting 2Abr is extracted as (e^ — 

e — )/2 for Abr < Ai, and (e++ - e — )/2 for Abr > Ai, 
at K. Figure [3ja) illustrates the zigzag potential mod- 
eling the transverse field. The calculated 2Abr versus 
the electric field E is shown in Fig. ^h). The depen- 
dence is linear with the slope of 10 /^eVnm/V. Since the 
field of 1 V/nm is produced by an electron charge 1 nm 
away from the graphene sheet, such fields are typical for 
graphene on a substrate. We expect that in realistic sit- 
uations the intrinsic and extrinsic spin-orbit couplings 
compete, making the topologies described in Fig.[2]hkely 
occurring in real samples. Since the extrinsic coupling 
depends linearly on the electric field, the topology is tun- 
able by gates. We also give the calculated magnitude of 
the graphene's dipole moment (the shift of the electron 
charge density): 0.0134 CA in unit cell. One may then 
relate the Bychkov-Rashba effect directly to the induced 
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dipolar moment; this should be particularly useful for es- 
timating the extrinsic splitting due to ad-atoms absorp- 
tion on graphene. 

Previous numerical estimates for the intrinsic split- 
ting 2Ai in graphene are rather controversial. The split- 
tingwas estimated to be in the range of 1 to 200 /leV 
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_ l26| . Kane and Mele [8| estimated the split- 

ting of 200 /LieV. This optimistic estimate was drastically 
reduced by Min et al. to the value of 1 ^eV, sup- 
ported by subsequent works [2J, |25[ . None of these stud- 
ies were fully first-principles. A density functional cal- 
culation of Boettger and Trickey [i^ , using a Gaussian- 
type orbital fitting function methodology, gave 50 //eV. 
Our result is about one half of that; the difference is likely 
due to the different ap prox imation schemes for spin-orbit 
coupling used in Ref. |26| and by us [2^ . Previous esti- 
mates for the extrinsic splitting, 2Abr, are 0.516 ^eV [§| 
and 133 /xeV [1^ per V/nm. No fully first-principles cal- 
culation of 2Abr, or the extrinsic effects in graphene in 
general, was reported thus far. 

What is the origin of the rather large, as compared to 
previous non-fully-first-principles results, intrinsic spin- 
orbit splitting in graphene? We calculate the spin-orbit 
coupling splitting of the 2p levels in carbon atoms to 
be 8.74 meV, using the Wien2k code. This splitting 
should be reflected in the splitting of the bands at the 
r point. Our calculation finds the splitting at the F 
point of 8.978 meV about 3 eV below the Fermi level, 
in close agreement with the atomic value. The bands at 
the K {K') points are formed mainly by pz orbitals whose 
magnetic quantum number is zero. The intrinsic splitting 
can be due to the coupling of the orbitals (forming the 
TT-bands) to either a bands or bands formed by higher or- 
bitals (d, /, ...). As argued already by Slonczewski [l^ 
using group theory, it is the d and higher orbitals that 
dominate the spin-orbit splitting at K{K'). A quali- 
tative argument for that was provided by McClurc and 
Yafet [l5| : Orbitals dxz and dyz can form Bloch states of 
the TT band symmetry at K{K'). Due to a finite overlap 
between the neighboring and dxz, dyz orbitals, the in- 
trinsic splitting is linearly proportional to the spin-orbit 
splitting of the d states (orbitals higher than d have a 
smaller overlap and contribute less). In contrast, due to 
the absence of the direct overlap between t he P z and a- 
band orbitals, the usually considered [1, [23, 24, 2^ spin- 
orbit splitting induced by the a-ir mixing depends only 
quadratically on the atomic spin-orbit splitting, giving a 
negligible contribution. 

In Fig. [3{c) we show the orbital-resolved densities of 
states. The Px,Py atomic character vanishes above —3 eV 
(the Fermi level ep is at zero). The density of states 
close to Ep comes predominantly from pz orbitals (tt- 
bands). Nevertheless, there is a finite contribution from 
d orbitals {dxz and dyz) which follows in shape that of 
Pz- The contributions from dx^-y^, dxy, and dz2 van- 
ish for energies above the — 3eV. To further confirm the 



symmetry arguments of Slonczewski, we have selectively 
removed orbitals from the calculation of the spin-orbit 
coupling contribution. Removing d and higher orbitals 
reduces the intrinsic splitting to 0.98 n eV, reproducing 
earlier non- first-principles calculations [2^, [3 ■ We 
conclude that the intrinsic splitting is dominated by d 
and higher orbitals, giving more that 96% of the splitting 
[30| . In contrast, we find that the extrinsic BR coupling 
is largely unaffected by the presence of d and higher or- 
bitals, demonstrating that this coupling is due to the ct-tt 
mixing. 

Finally, we calculate the spin-orbit splitting for what 
we call a 'mini-ripplc' configuration (Fig.[3|), in which ev- 
ery other atom is displaced by an amount A transverse 
to the sheet. This calculation should give an indication 
of what to expect for larger-scale ripples that occur in 
graphene on a substrate 27 1 or free standing [gH; such 
large scales are out of the scope for our methods. The 
mini-ripple exhibits a gap that grows quadratically with 
increasing A, as seen from Fig. 2) Removing d and higher 
orbitals from the calculation of the spin-orbit splitting, 
the initial gap reduces to about 1 ^ eV but the overall 
growth remains largely unchanged: the rippling-induced 
gap is almost solely due to ct-tt mixing. Furthermore, 
the displacements of less than 1% have no significant ef- 
fects on the intrinsic spin-orbit splitting. For an effective 
description of the gap opening in the mini-ripple we con- 
sider a more general extrinsic case. The contribution of 
a transverse electric field is twofold. First, there is a 
pseudospin splitting, as the two sublattices are no longer 
equivalent (the two triangular sublattices have a different 
potential). Second, the Bychkov-Rashba effect appears. 
We find that the following formula describes the resulting 



spectrum: e 



K 

^1/ 



(2Abr)= 



where Aa is the electric field induced gap due to the rip- 
pling (Aa oc E). 

In summary, we have shown that for realistic elec- 
tric fields the graphene band structure exhibits remark- 
able tunable topologies. Our first-principles calculation 
gives strong support for the effective spin-orbit coupling 
Hamiltonian models, making them highly accurate ana- 
lytical tools to investigate the physics of graphene. 
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FIG. 1: Graphene's essentials. Bottom-up: Carbon atoms form a honeycomb lattice with two atoms in the unit cell. The first 
Brillouin zone of the reciprocal lattice contains two nonequivalent Dirac points, K and K' . The relevant states at the Fermi 
level form two touching cones with the tips at K{K'). 




FIG. 2: Band structure topologies of graphene. Transverse electric field drastically changes the topology of the bands near 
K{K'). (a) Zero electric field, (b) Electric field of magnitude E = l.OV/nm. (c) E = 2.44 V/nm. (d) E = 4.0 V/nm. The 
first-principles results are represented by circles (the Fermi level is at zero). The curves are fits to the analytical model (see 
the main text). The spin branch = f is shown in red, /i = — f in blue. The calculated Fermi velocity is up = 0.833 X iO^ m/s. 
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FIG. 3: Bychkov-Rashba-type splitting in graphene and orbital-resolved density of states, (a) Transverse electric field is 
modeled with a zigzag potential, (b) Bychkov-Rashba spin-orbit induced splitting at K{K') as a function of the electric field. 
The slope is 9.9 neV per V/nm. (c) Projected density of states to particular atomic orbitals. The Fermi level is at zero. The 
d-character is enhanced by a factor of ten. 
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FIG. 4: Calculated spectral gap at K{K') in the mini-ripple configuration (inset), as a function of the relative corrugation 
strain with respect to the lattice constant a. Neglecting d and higher orbitals results in an almost constant shift, proving that 
the main effects come from a-n hybridization. 



